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A  scanning  acoustic  microscope  that  uses  a  novel  lens  having  a  slot-shaped 
aperture  is  used  to  investigate  the  quantitative  measurements  of  properties  of  isotropic 
and  anisotropic  materials  and  also  the  characterization  of  surface-  breaking  cracks.  The 
slot  lens  has  a  point  focus  and  yet  launches  leaky  Rayleigh  waves  in  preferential 
directions.  This  allows  surface  wave  velocity  measurements  of  anisotropic  materials  as 
well  as  interrogation  of  cracks  at  nearly  normal  incidence  for  determination  of  the  crack 
reflection  coefficient.  Such  a  lens  can  also  be  used  for  imaging  of  polycrystalline 
material  surfaces  where  contrast  depends  solely  on  Rayleigh-wave  anisotropy. 
Diffraction-corrected  geometric  theory  is  employed  to  model  the  response  of  both  the 
slot  and  circular  aperture  lenses  to  cracked  as  well  as  defect-free  sample  surfaces. 
Experimental  measurements  were  made  to  test  the  theoretical  predictions. 
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Chapter  I 
Introduction 


1.1  Acoustic  microscope 

The  acoustic  microscope  used  in  this  work  consists  of  a  fused  quartz  lens  rod 
which  has  a  transducer  bonded  to  one  end  and  a  concave  spherical  lens  ground  in  the 
other  end  (Fig.  1.1).  To  make  a  measurement  at  a  particular  location,  a  tone-burst  signal 
of  49  MHz  and  0.5  jisec  duration  is  applied  to  a  piezoelectric  transducer.  The  transducer 
excites  a  compressional  acoustic  wave  which  travels  down  the  lens  rod.  The  curved 
surface  of  the  lens  causes  focusing  in  the  coupling  fluid  (usually  water)  at  or  near  the 
surface  of  the  sample.  The  resulting  incident  wave  travels  as  a  conical  beam,  thus 
illuminating  the  sample  with  rays  that  cover  a  wide  range  of  angles.  The  tone-burst 
signal  reflected  from  the  surface  of  the  sample,  which  now  contains  information  about 
the  sample,  is  refracted  at  the  lens  surface  into  the  lens  rod  and  ultimately  excites  the 
same  transducer.  Time-gating  is  employed  in  the  receiver  circuit  to  isolate  the  desired 
echo  signal  and  exclude  internal  rod  reverberations  and  other  undesired  signals  that  may 
be  present.  The  lens  can  be  mechanically  scanned  along  any  of  the  three  axis,  thus 
allowing  various  types  of  measurements,  including  imaging,  to  be  performed. 

It  has  been  shown  by  Weglein  [1],  Bertoni  [2],  and  others  that  when  the  focal 
point  of  the  lens  lies  below  the  sample  surface,  the  signal  returning  from  the  sample  can 
be  viewed  as  consisting  of  two  contributions:  a  geometrically  scattered  part  and  a  leaky 
wave  part.  The  geometric  contribution  results  from  specular  reflection  from  the  sample 
surface.  The  leaky  contribution  results  from  those  incident  rays  that  strike  the  sample  at 
the  Rayleigh  critical  angle,  thus  launching  Rayleigh  surface  waves  which  propagate 
along  the  sample  surface  and  radiate  (leak)  energy  back  into  the  water  in  the  form  of 
P-waves  (  Fig.  1.2).  Some  of  these  leaked  rays  are  refracted  by  the  lens  toward  the 
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Transducer 


Fig.  1.1  Reflection  mode  acoustic  microscope  with  sample  at  focus. 
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Fig.  1.2  Ray  diagram  of  the  acoustic  lens,  defocused  toward  the  sample,  showing  the 
specularly  reflected  axial  ray  and  the  excited  and  re-radiated  leaky  surface 
rays. 
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transducer  where  they  become  part  of  the  total  transducer  voltage. 


1.2  The  V(z)  measurement 

An  important  type  of  measurement,  called  the  acoustic  material  signature  or  V(z), 
involves  the  variation  in  amplitude  of  the  signal  reflected  from  the  sample  surface  as  the 
lens  is  defocused  toward  the  sample  along  the  vertical  or  z-axis.  The  defocus  distance  z 
is  measured  from  the  front  focal  plane  of  the  lens.  The  rays  comprising  the  geometric 
and  leaky  contributions  travel  along  different  paths  which  results  in  the  two  components 
arriving  at  the  transducer  with  different  phases.  This  phase  difference  is  a  function  of  z, 
which  leads  to  an  interference  pattern  in  V(z)  resulting  from  the  geometric  and  leaky 
contributions  gradually  going  in  and  out  of  phase  as  the  lens-sample  separation  is  varied. 
A  typical  V(z)  measurement  made  using  as  a  sample  a  glass  microscope  slide  is  shown  in 
Fig.  1.3.  The  separation  between  local  minima  is  directly  related  to  the  Rayleigh  wave 
velocity,  which  can  thus  be  extracted  from  the  V(z).  This  velocity  is  strongly  dependent 
upon  the  elastic  properties  of  the  sample  and  thus  the  V(z)  measurement  becomes  a 
valuable  tool  in  the  evaluation  of  material  properties  . 

1.3  The  V(x)  measurement 

Another  important  type  of  measurement,  called  V(x),  records  the  transducer 
voltage  as  the  lens  is  scanned  laterally  at  constant  defocus  height  z  above  the  sample. 
When  the  sample  surface  contains  a  discontinuity,  such  as  a  thin  crack,  the  transducer 
voltage  is  altered  since  the  defect  may  interfere  with  the  surface  wave  propagation.  The 
usual  leaky  contribution  is  reduced  and  another  contribution  to  the  total  detected  signal 
arises  due  to  the  surface  waves  that  are  scattered  from  the  crack.  Some  of  the  rays 
reflected  from  the  crack  or  transmitted  through  it  are  shown  in  Fig.  1.4.  The  strength  of 


4 


V(z) 

Relative 

Units 


5 


Fig.  1.3  Measured  transducer  output  voltage  V  as  a  function  of  defocus  distance  z  for  a 
glass  sample  and  a  circular  aperture  lens 


crack 


I 


Fig.  1.4  Rayleigh  wave  interaction  with  a  crack.  A  surface  wave  normally  incident  on 
the  crack  from  the  left  results  in  reflected  and  transmitted  surface  waves 
propagating  away  from  the  crack. 
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Fig.  1.5  Measured  transducer  output  voltage  V  as  a  function  of  lateral  displacement  x 
for  a  glass  slide  containing  a  deep  surface  breaking  crack  at  x  =  0. 


the  scattered  signal  is  directly  proportional  to  the  crack  surface  reflection  coefficient. 
The  propagation  path  length  of  the  component  scattered  from  the  crack  varies  with  crack 
displacement  x.  This  results  in  interference  between  this  component  and  the  geometric 
contribution,  whose  phase  remains  constant  for  all  x.  A  typical  V(x)  measurement  of  a 
cracked  glass  slide  is  displayed  in  Fig.  1.5.  Angel  and  Achenbach  [3]  have  shown  that 
the  surface  wave  reflection  coefficient  R<.  at  the  crack  is  related  to  crack  depth.  Thus  if 
Rc  is  determined  from  V(x),  it  may  be  possible  to  determine  the  crack  depth,  which  is  of 
interest  to  non  destructive  evaluation  (NDE). 

1.4  The  slot  lens 

Conventionally,  a  point-focus  lens  with  a  circular  aperture  is  used  for  V(z) 
measurements.  When  this  type  of  lens  is  defocused  toward  a  sample,  surface  waves  are 
launched  in  all  directions  on  the  sample  surface  which  makes  this  lens  unsuitable  for 
meaningful  measurements  of  anisotropic  materials  as  the  Rayleigh  wave  velocity  varies 
with  the  direction  of  propagation.  In  addition,  a  circular  lens  cannot  readily  be  used  to 
measure  surface  wave  reflection  coefficient  of  cracks,  since  this  coefficient  is  also 
strongly  dependent  on  the  angle  of  incidence  of  the  surface  wave  with  the  crack. 

A  line-focus  (cylindrical)  lens  has  been  used  by  Kushibiki,  et.  al.  [4]  with 
anisotropic  materials  to  launch  surface  waves  along  a  particular  direction  of  the  sample 
surface.  However,  since  it  focuses  to  a  line  and  not  a  point,  it  has  good  resolution  in  only 
one  direction  and  is  thus  unsuitable  for  imaging  polycrystalline  materials.  This  lens  can 
also  be  used  for  determination  of  the  crack  reflection  coefficient,  however,  again,  its  poor 
lateral  resolution  requires  a  very  long  straight  uniform  crack  which  may  not  be  feasible 
experimentally  or  encountered  in  practice. 

To  conduct  measurements  of  anisotropic  crystalline  materials,  graphite-epoxy 
composites,  as  well  as  defect  characterization,  a  new  kind  of  lens  was  proposed  [5].  In  a 
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Bottom  View 


Fig.  1.6  Slot  aperture  acoustic  lens 
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novel  approach,  part  of  the  aperture  of  the  new  lens  was  blocked  by  an  acoustically 
absorbent  material  so  as  to  leave  only  a  slot-shaped  aperture  exposed.  This  lens,  called 
the  slot  lens  (Fig.  1.6),  can  launch  and  receive  leaky  surface  waves  only  along  a  narrow 
range  of  directions,  thus  allowing  anisotropy  measurements  to  be  performed,  with  only 
somewhat  degraded  resolution.  It  can  also  be  used  to  interrogate  short  cracks  under 
nearly  normally  incident  conditions  which  allows  the  determination  of  the  crack 
reflection  coefficient,  and  thus,  possibly,  the  crack  depth. 

To  assess  the  lens  performance  it  is  necessary  to  have  a  theoretical  model  to  be 
able  to  predict  the  lens  response  to  both  cracked  and  defect-free  surfaces.  Atalar  [6]  used 
the  angular  spectrum  method  to  develop  a  theory  for  predicting  V(z).  It  takes  into 
account  field  diffraction  both  in  the  lens  rod  and  in  water  and  treats  the  concave  surface 
of  the  lens  as  a  thin  lens.  Bertoni  [2]  has  proposed  an  alternate  approach  that  uses 
geometrical  ray  techniques  to  model  the  microscope.  While  the  latter  theory  ignores  the 
effects  of  diffraction  which  produces  poor  agreement  near  focus,  it  provides  a  more 
physical  visualization  as  well  as  computationally  simpler  approach. 

It  is  also  important  to  model  the  V(x)  for  the  purpose  of  extracting  the  crack 
surface  wave  reflection  coefficient.  Somekh,  et.  al.  [  7]  extended  Atalar's  [6]  approach 
to  cracked  surfaces  measured  with  a  cylindrical  lens.  Ahn,  et.  al.  [8]  have  used  a  Green's 
function  solution  implemented  by  numerical  boundary  element  method.  None  of  these 
approaches  involve  the  reflection  coefficient  of  the  crack  in  a  simple  way  that  permits  a 
direct  determination  of  this  coefficient. 

In  this  work  geometrical  ray  techniques  are  extended  to  include  diffraction  effects 
in  the  lens  rod  and  are  used  to  model  the  V(z)  and  the  V(x)  for  both  circular  and  slot 
aperture  lenses.  A  simple  technique  is  developed  to  extract  the  crack  reflection 
coefficient  from  the  V(x). 
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Chapter  H 

Theoretical  modelling  of  the  V(z)  for  the  circular  aperture  lens 


A  review  of  two  theoretical  models  that  predict  V(z)  for  a  circular  aperture  lens 
are  presented:  the  angular  spectrum  approach  of  Atalar  [6]  and  the  ray  method  of  Bertoni 
[2].  A  modification  of  the  ray  method  is  then  proposed  that  accounts  for  diffraction 
effects  in  the  lens  rod.  Finally,  a  method  of  extracting  the  Rayleigh  wave  velocity  from 
the  V(z)  is  described. 

2.1  The  angular  spectrum  approach 

The  following  is  a  brief  overview  of  the  angular  spectrum  model  of  the 

microscope  V(z)  response,  as  developed  by  Atalar  [6]. 

The  transducer  radiates  an  acoustic  field  u0+(x,y)  that  propagates  down  the  lens 
rod,  undergoes  diffraction  and  is  defined  at  the  back  focal  plane  of  the  lens  as  the  field 
ui+(x,y).  A  crucial  assumption  is  that  the  spherical  lens  is  assumed  to  be  thin  and,  thus, 
the  field  U2+(x,y)  at  the  front  focal  plane  may  be  expressed  in  terms  of  the  Fourier 
transform  of  U!+(x,y).  The  phase  variation  produced  in  propagating  from  the  front  focal 
plane  to  the  plane  of  the  sample,  located  a  distance  z  from  the  front  focal  plane  is 
included  by  multiplying  the  angular  spectrum  of  the  field  U2+(x,y)  by  the  phase  factor 
exp  (i  kz  z).  The  resultant  angular  spectrum  is  multiplied  by  the  water-sample  reflection 
coefficient  #(kx,  ky)  and  multiplied  again  by  the  above  phase  factor,  which  results  in  the 
reflected  spectrum  U2'(kx,ky)  at  the  front  focal  plane.  The  field  U2‘(x,y)  is  then  Fourier 
transformed  again  to  obtain  the  reflected  field  uj+(p)  at  the  back  focal  plane.  The 
transducer  response  voltage  V  is  found  by  integrating  the  product  of  the  acoustic  fields 
defined  at  the  back  focal  plane,  i.e. 
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(2.1.1) 


oo  oo 

V(z)  =  /  J  U!+(x,y)  uf (x,y)  dx  dy 


Substituting  uf(x,y)  found  by  carrying  out  the  above  steps  leads  to 


oo  oo 

V(z)  =  f  f  Uf(x,y)  Uf( -x, -y)  ^(kx,  ky)  PiCx.y)  P2(x,y)  exp(i2kz  z)  dx  dy  (2.1.2) 

-oo  -oo 


In  evaluating  the  above  integral,  kx,  ky  and  kz  are  expressed  in  terms  of  x  and  y  by  using 
the  following  relations 

kx  =  kw  sin  9  cos  <]> 
ky  =  kw  sin0sin<}> 
kz  =  kw  cos  9 


0  =  tan  P  /  (f  -  R  +  “\/r2  -  p2) 
(J)  =  tan y/x 


In  these  relations  kw  and  f  are  the  wavenumber  and  focal  length,  respectively,  in  water, 
and  R  is  the  lens  radius  of  curvature.  The  pupil  functions  Pj(x,y)  and  P2(x,y)  account  for 
two-way  propagation  through  the  lens  aperture,  and  include  both  the  finite  extent  of  the 
aperture  as  well  as  the  lens-water  and  water-lens  transmission  coefficients  tLW  and 
Twl>  respectively.  In  the  usual  Fourier  optics  approximation  [9],  these  functions  are 
evaluated  at  the  lens  aperture  and  projected  onto  the  back  focal  plane  in  unmodified 
form,  thus  neglecting  diffraction  effects  on  the  pupil  functions  between  these  two  planes. 

The  expression  (2.1.2)  for  V(z)  can  be  evaluated  numerically.  The  results  of  such 
calculations  are  plotted  against  z  in  Fig.  2.1  and  2.2  for  teflon  and  glass  samples, 
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the  angular  spectrum  approach  (dotted). 
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the  angular  spectrum  approach  (dotted). 


respectively.  The  water-teflon  interface  does  not  support  surface  waves,  while  the  water- 
glass  interface  does,  as  is  evidenced  by  the  interference-like  structure  in  its  V(z). 
Comparisons  are  made  with  the  measured  values  for  these  materials  in  Fig.  2. 1  and  2.2  . 

Strictly,  T^l  is  also  a  function  of  sample  defocus  z,  as  the  transmission 
coefficient  depends  on  the  angular  spectrum  of  the  reflected  field.  In  Fig.  2.3  a  V(z)  for 
teflon  calculated  using  T^Cz)  is  compared  with  a  V(z)  calculated  while  assuming 
tWl(z)  =  Twl(z=0)-  As  the  z  variation  of  apparently  does  not  significantly  effect 
the  V(z),  Twl(z=0)  will  be  used  for  all  z,  This  simplification  was  used  in  producing  Fig. 
2.1  and  2.2.  Both  TLW  and  T^  are  calculated  using  the  method  of  acoustic  equivalent 
networks  for  planar  interfaces  by  Oliner,  et.  al.  [10]. 

While  the  angular  spectrum  approach  gives  impressively  good  agreement  with  the 
experimental  results,  it  is  difficult  to  gain  intuitive  insight  into  the  physical  processes 
involved  in  determination  of  V(z)  from  the  integral  (2.1.2). 
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Fig.  2.3  V(z)  of  Teflon  calculated  using  (a)  T^Cz)  (solid)  and 
(b)  assuming  T^(z)  =  T^Cz  =  0). 


2.2  The  Ray  approach 


The  following  is  a  brief  overview  of  the  ray  model  of  the  V(z),  as  developed  by 
Bertoni  [2]. 

The  methods  of  ray  optics  used  here  are  well  described  by  James[ll].  To 
conserve  power,  the  field  amplitude  for  spherically  convergent  or  divergent  beams  varies 
inversely  as  the  square  root  of  the  cross-sectional  area  of  the  beam.  For  cylindrically 
focussed  beams  this  variation  becomes  inversely  proportional  to  the  square  root  of  the 
distance  from  the  focus.  Furthermore,  while  propagating  through  focus,  the  beam  suffers 
a  phase  change  of  -7t  for  spherical  (point)  focus,  and  -tc/2  for  a  cylindrical  focus. 

The  transducer  response  voltage  V(z)  consists  of  two  parts:  the  specular 
(geometric)  contribution  Vq(z)  and  the  leaky  contribution  Vl(z) 

V(z)  =  VG(z)  +  VL(z).  (2.2.1] 

These  two  components  are  modelled  separately: 

2.2.1  The  geometric  component  Vq 

Neglecting  diffraction,  the  transducer  illumination  of  the  back  of  the  lens  is  a 
uniform  compressional  acoustic  field.  Assuming  axial  symmetry,  the  field  returning  from 
the  sample  and  striking  the  plane  of  the  transducer  is  found  to  be: 

Iful  r  ikw  "I 

■Tlw  T\vl  ^(0)  exp(i2kw(D/n+f+z)  |f+2zl  ID-ul  exP  L  2n(D-u)  J  ’  (2.2.2 ) 

where  u  =  f(f+2z)/2nz,  D  is  the  distance  between  the  transducer  and  the  lens  aperture,  n 
is  lens  material  index  of  refraction  relative  to  water  (n  =  vj/vw),  and  f  is  the  lens  focal 
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length  in  water,  which  can  be  expressed  in  terms  of  the  lens  radius  of  curvature  R 


(2.2.3) 


The  lens  transmission  coefficients  Tl^t  and  T^l  as  well  as  the  sample  reflection 
coefficient  IS  are  assumed  here  to  be  constant  for  all  angles  of  incidence.  The  derivation 
of  (2.2.2)  will  be  considered  in  more  detail  in  the  next  section.  The  voltage  Vq  is 
obtained  by  integrating  the  field  (2.2.2)  over  the  transducer  and  normalizing  by  the 
transducer  area,  resulting  in 


_  2k  n(D-u) 

VG  =  -Tlw  Twl  0(0)  exp  [i  2  kw  (D/n+f+z)]  jj- 

[  i  kw  rT2  1 

•  LexP5^)  ■  U  1 


(2.2.4) 


where 

rT  =  ^a»  zmin >z  or  z>  zmax 

=  Ra(l  -  2z  (nD  -  f)  /  f),  Zujin  <  z  <  0 

=  ®-a  ®-ul  /  u,  0  <  z  <  Zjjjgx 

£  (R^Ra  ~  D 

zmin  -  _  2n(D-f/n) 

f2(l+RT/Ra) 

zmax  -  2n[D-(l+R’p/Ra)f/n  * 


(2.2.5) 


[2.2.6] 


Here  Ra  is  the  aperture  radius.  The  expression  (2.2.4)  was  evaluated  numerically  for  a 
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the  geometric  ray  approach  (dotted). 


teflon  sample  and  plotted  in  Fig.  2.4  along  with  the  measured  V(z)  of  teflon.  Since 
Rayleigh  waves  cannot  be  excited  on  teflon  from  the  water,  there  is  no  leaky 
contribution.  One  may  observe  in  Fig.  2.4  that  the  agreement  with  the  experiment 
deteriorates  as  Izl  approaches  focus  as  is  expected  when  diffraction  is  neglected. 


2.2 2  The  leaky  component  VL(z) 


First,  the  behavior  of  VL(z)  for  distances  closer  than  geometric  focus  (z  <  0)  is 
considered.  Leaky  waves  are  excited  by  rays  that  are  incident  on  the  sample  surface  at 
angle  0r  or,  equivalently,  strike  the  lens  surface  at  p  =  Pr  (derived  in  Appendix  A), 
from  the  lens  axis,  as  shown  in  Fig.  2.5.  The  resulting  surface  wave  propagates  along  the 
surface  leaking  a  parallel  family  of  rays  at  0R  with  respect  to  the  normal  to  the  surface. 
Since  the  lens  is  axially  symmetric,  these  leaky  waves  are  launched  in  all  directions  on 
the  sample  surface  passing  through  the  center  of  symmetry  defined  by  the  lens  axis.  This 
center  of  symmetry  thus  becomes  a  cylindrical  focus  for  the  surface  waves.  In  any  radial 
plane,  one  of  the  leaked  rays,  called  the  principal  ray,  radiates  in  such  a  way  as  to  cross 
the  lens  aperture  plane  at  p  =  p-p  and  travel  to  the  transducer  along  a  path  that  is  the 
mirror  image  of  the  incident  ray.  A  bundle  of  rays  neighboring  the  principal  ray  focus  to 
a  point  in  the  radial  plane  located  in  the  lens  rod  a  distance  fx  (derived  in  Appendix  A) 
above  the  aperture.  Due  to  the  symmetry  of  the  system  with  respect  to  the  lens  axis  (z) 
the  entire  family  of  leaky  rays  forms  a  ring  focus  of  radius  Pr  between  the  lens  surface 
and  the  transducer.  The  response  voltage  VL  is  obtained  by  integrating  the  resulting  field 
over  the  transducer  and  normalizing  to  the  transducer  area,  resulting  in 


[ 


2<*L  K  (V*  f3  si"  9.)1'2  ]  : 
- Wiw  Je  ’K 


VL=  -  Tlw  Twl  exp(i2kw(D/n+f)) 

•  exp(i2kw(D/n+f+lzlcos0R)  •  exp(-2cxjzl  tan  0r), 


(2.2.7) 
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where 


Nt 

K  =  2e^  /  [exp  i(N-NR)2+  exp(-ijt/2)  exp  i(N+NR)2  ]  a/n  dN  (2.2.8) 

^  0 


In  (2.2.8), 


Nt  = 


RT  and  Nr  = 


PR- 


The  expressions  (2.2.7)  and  (2.2.8)  can  be  evaluated  numerically  to  yield  Vl(z)  for  z  <  0. 
For  z  >  0  the  leaky  rays  miss  the  transducer  and,  therefore: 


VL(z)  =  0  for  z  >  0 


(2.2.9) 


2.2 3  The  total  voltage  V(z) 

The  total  voltage  V(z)  is  found  by  adding  Vq  and  Vl«  The  result  of  such 
calculation  for  a  glass  sample  is  compared  in  Fig.  2.6  with  the  measured  V(z).  It  is  seen 
that,  as  in  the  case  of  teflon  (Fig.  2.4),  the  agreement  with  the  experiment  deteriorates  as 
Izl  approaches  focus. 
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Fig  2.5  Cross-sectional  view  of  the  lens  showing  a  bundle  of  leaky  rays  forming  a  focus 
in  the  lens  rod. 
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the  geometric  ray  approach  (dotted). 


2.3  Diffraction  correction  to  the  ray  model 


2.3.1  Introduction 

The  theoretical  V(z)  obtained  with  the  angular  spectrum  approach  agrees  much 
better  with  the  experimental  data  then  the  ray  model  V(z),  especially  near  focus  (  Izl 
small),  as  can  be  seen  in  Figs.  2.4  and  2.6.  Since  the  angular  spectrum  model  takes 
diffraction  into  consideration,  it  may  be  suspected  that  diffraction  effects  in  the  lens  rod 
play  an  important  role  in  forming  V(z).  Therefore  it  is  desirable  to  compensate  for  these 
effects  in  the  ray  model  of  V(z). 

Assuming  the  transducer  radiates  as  a  uniform  piston  source,  the  diffracted  field 
U  at  the  back  of  the  lens  may  be  found  by  the  method  of  Tjotta  and  Tjotta  [12].  The 
amplitude  and  the  phase  of  this  field  are  plotted  as  a  function  of  p  in  Figs.  2.7a  and  2.7b. 

Transmission  and  refraction  at  the  lens-water  interface  will  now  be  modelled  by 
geometrical  ray  theory  which  assumes  that  the  field  behaves  locally  as  a  plane  wave  and 
treats  the  phase  and  amplitude  variations  of  the  field  separately.  The  gradient  of  the 
phase  of  the  field  at  a  point  determines  the  direction  of  the  local  wavevector  and,  thus, 
the  local  angle  of  incidence.  The  direction  of  the  refracted  ray  may  then  be  found  using 
Snell's  Law,  thereby  defining  the  location  of  focus.  The  angle  of  incidence  also 
influences  the  lens-water  transmission  coefficient  tLW>  thus  effecting  the  amplitude 
profile  of  the  transmitted  field. 
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2.3.2  Approximation  of  phase  variation  of  field  U(p) 


The  local  wavevector  k  may  be  expressed  as 

K(£)  =  V  <)>(r),  (2.3.1) 

where  <)>  is  the  spatial  component  of  the  phase  factor,  i.e. 

U(e)  =  IU(r)l  eW  ,  (2.3.2) 

In  the  lens  rod  which  has  axial  symmetry  about  the  z  axis 


E  =  p£!o  +  zz0, 


(2.3.3) 


and  l]t(l)l  =  kL  =  co  /vl  is  the  wavenumber  in  the  lens  medium.  From  (2.3.1)  the 
component  of  along  the  p  direction  is 


(2.3.4) 


The  angle  of  propagation  0,  measured  with  respect  to  the  z-  axis  is  then  found  as 


0  =  sin*1  (kp  /  kL)? 


(2.3.5) 


The  angle  0,  found  from  the  phase  of  the  field  by  numerically  carrying  out  the 
above  procedure  is  plotted  as  a  function  of  p  in  Fig.  2.8.  For  lens  considered  here  the 
maximum  angular  deviation  of  the  incident  ray  from  the  vertical  axis  is  about  1.3°.  It  is 
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therefore  reasonable  to  assume  that  the  acoustic  field  incident  on  the  lens  surface  has  no 
cross-  sectional  phase  variation  and  thus  strikes  this  surface  as  a  well-  collimated  beam. 
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Fig  2.7b  Relative  phase  profile  of  field  U(p)  illuminating  the  back  focal  plane. 
Calculated  using  the  theory  of  Tjotta  and  Tjotta  [12]. 
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vertical  of  rays  illuminating  the  b 


2.3 3  Approximation  of  amplitude  variation  of  field  U(p) 


The  diffracted  field  U(p)  is  transmitted  through  the  lens  surface  into  water  with 
a  transmission  coefficient  TLW(p).  The  field  reflected  from  the  sample  surface  should 
similarly  be  scaled  by  the  water-lens  reflection  coefficient  TyyL(P)-  As  has  been  justified 
in  the  angular  spectrum  approach,  T^l  is  assumed  to  be  independent  of  z  and  the  field 
U(p)  TLW(p)  may  therefore  be  multiplied  by  T^l  prior  to  scattering  from  the  sample 
surface.  Normalizing  by  T^fO)  TLW(0),  the  resulting  field  variation  Ueff(p)  is  given 
by 


Ueff(p)  = 


U(p)  TL\y(p)  T\vl(P)  exp(kw  D/n) 
Twl(O)  Tlw(^) 


(2.3.6) 


The  amplitude  of  Uef{(p)  is  plotted  in  Fig.  2.9.  The  aperture  pupil  function  Pi  is  defined 
as 


Pl(p)  =  circ  (p) .  (2.3.7) 

It  is  possible  to  approximate  the  amplitude  of  Ueff  with  a  superposition  of 
Gaussians 


Ueff(p)=  exp(ikw  D/n)  E;  A.  exp(^) 


(2.3.8) 


The  result  of  fitting  with  two  Gaussians: 


UefiCp)  =  exp(ikw  D/n)  £,•  Af  exp(^)  i=l,2 


(2.3.9) 
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is  plotted  in  Fig.  2.9  for  comparison  with  the  exact  field  Ueff.  The  amplitudes  A,-  and 
widths  B,-  of  these  Gaussians  were  found  by  trial  and  error  to  be  =  1.35,  Bj  =  1000 
jxm,  A2  =  .58,  and  B2  =  2700  pm  for  the  lens  considered  here  (see  Appendix  B  for  lens 
dimensions).  It  is  possible  to  use  a  least  squares  fit,  as  well  as  more  terms  in  expansion 
(2.3.8)  to  improve  the  approximation  to  IUeffl.  However,  it  appears  from  Fig.  2.9  that  for 
the  lens  considered  here  the  use  of  these  two  Gaussians  is  a  sufficiently  good 
approximation  of  lUeff)  and,  therefore,  this  representation  will  be  used  here. 


2.3.4  Calculation  of  the  transducer  response  Vq(z) 

It  is  first  important  to  summarize  the  microscope  ray  analysis  as  developed  by 
Bertoni  [2]  for  the  case  of  uniform  illumination  of  the  lens.  The  source  field  U  is 
transmitted  through  the  lens  and  focused  to  a  point  at  the  front  focal  plane.  The  sample 
surface  acts  as  a  plane  mirror  located  a  distance  z  away  from  the  front  focal  plane,  and, 
therefore,  creates  an  image  of  the  focal  point  at  a  distance  f+2z  below  the  lens  aperture 
(Fig.  2.10).  The  field  suffers  a  phase  change  of  7t  due  to  passage  through  focus.  The 
resulting  field,  which  now  travels  away  from  the  sample,  is  imaged  by  the  lens  to  form 
an  effective  point  source  illuminating  the  transducer.  The  location  u  of  this  source  can 
be  found  from  the  lens  imaging  formula  as 


u  = 


f(f+2z) 

2nz 


(2.3.10) 


The  sample  reflection  coefficient  is  assumed  to  be  constant,  i.e.  $(p)  =  $(0),  which  is 
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Ueff  (p)| 
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Fig.  2.9  Amplitude  profile  of  effective  field  Ueff(p)  illuminating  the  back  focal  plane: 

(a)  Calculated  using  the  theory  of  Tjotta  and  Tjotta  [12]  (solid)  and 

(b)  two-Gaussian  fit  (dashed). 


h-  Ri  H 


Fig.  2.10  Ray  structure  for  specularly  reflected  component  VG  for  z  <  0. 
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strictly  true  at  normal  incidence.  Therefore  the  scattered  field  Ur(p)  illuminating  the  back 
of  the  lens  is  given  by: 


Tr\v(0)  T'y^tl(O)  m exp  (2i kyy (f+z)) 
Ur(P)= - Ml - 


(2.3.11) 


where  Mj  is  the  magnification  /  reduction  scaling  factor  due  to  ray  spreading  in  the 
water: 


M1  = 


If  +  2zl 
f 


The  aperture  pupil  functions  Pi  and  P2  are  defined  as 


?l(p)  =  p2(P)  =  circ  (P^a) 

circ(x)  =  1,  0  <  x  <  1 

=  0 ,  x  >  1. 


(2.3.12) 


(2.3.13) 


The  transducer  response  VG  can  be  found  from  (2.1.1).  Here  the  effective  source  field 
ui+(p)  is  Ueff(p)  and  the  scattered  field  uf(p)  is  Ur(p),  given  by  (2.3.11).  Using  the 
Gaussian  expansion  (2.3.9)  for  Ueff(p)  and  normalizing  by  the  transducer  area: 
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vG=- 


tLW(0)  tWL(0)  ®(0)  exP  (2i  kw  (D/n+f+z)) 


M^R2 


2?c 


Rn 


•  /  “P  (m^uT)  ^  A< A,  «p(^)  exp((B;M,)=)  P'(m7)  P‘<P)  P  dp 


0 


'•Mj' 


(2.3.14) 


The  product  of  sums  in  (2.3.14)  can  be  expanded: 


Ai,-A.A.  ,  (Bg*  "  Bj2  +  (B.  Mx)2 


(2.3.15) 


where  q  ranges  from  1  to  4  (assuming  two  Gaussians  are  used  in  (2.3.8)). 


Substituting  (2.3.15)  into  (2.3.14)  and  carrying  out  the  integration  gives 


vG  =  - 


tLw(0)  twl(Q)  m 
Rt2  Mj 


exp  (2i  kw  (D/n+f+z)) 


nBlft-u) 
i  kw  Bl,-2  -  2n(-u) 


,(2.3.16) 


where 


rT  =  Ra 


=  R 


Mj  <  1 

Mx  >  1 
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As  can  be  seen  from  the  above,  the  advantage  of  Gaussian  expansion  (2.3.8)  of 
IUeff(p)l  is  that  this  representation  allows  the  integration  to  be  carried  out  analytically. 

The  diffraction-corrected  voltage  Vq(z),  calculated  by  using  (2.3.16),  is  plotted 
versus  z  in  Fig.  2.11  along  with  Vq(z)  obtained  by  the  angular  spectrum  approach  and  in 
Fig.  2.12  for  comparison  with  the  experimental  Vg(z)  obtained  for  a  Teflon  sample, 
taking  advantage, of  the  fact  that  Teflon- water  interface  does  not  support  leaky  modes 
and  thus  only  the  specular  component  of  V(z)  exists.  In  Fig.  2.13  three  V(z)  curves  are 
plotted:  Vq(z)  calculated  assuming  uniform  lens  illumination  by  the  ray  method,  VG(z) 
corrected  for  diffraction  and  the  actual  Teflon  data.  Note  the  significant  improvement  in 
agreement  between  the  experimental  and  theoretical  curves  obtained  by  including 
diffraction  effects  in  the  lens  rod. 

2.3.5  Calculation  of  the  leaky  contribution  Vl(z) 

Now  the  effect  of  diffraction  in  the  lens  rod  on  the  leaky  contribution  VL(z)  will 
be  considered.  As  discussed  in  previous  section,  VL  arises  from  the  leaky  surface  waves 
excited  on  the  sample  surface  by  those  rays  from  the  transducer  that  cross  the  lens 
surface  at  p  =  pR.  Thus  the  amplitude  of  Vl  is  directly  proportional  to  the  amplitude  of 
the  diffracted  source  field  at  p  =  pR,  i.e.  IU(p)l.  Since  VL  in  (2.2.7)  was  calculated  using 
IU(p)l  =  1  for  all  p,  it  would  therefore  be  proper  to  scale  Vl  by  IU(pR)l  thereby  taking 
diffraction  in  the  lens  rod  into  account. 

The  amplitude  of  resulting  Vl(z)  is  shown  in  Fig.  2.14.  If  desired,  the  behavior  of 
Vl(z)  near  z  =  0  can  be  modelled  better  by  implementing  a  correction  described  by 
Bertoni  in  [2]. 
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(b)  the  diffraction-corrected  ray  approach  (dotted). 


V(z) 

Relative 

units 
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using  the  diffraction-corrected  ray  approach  (dotted). 


V(z) 

Relative 

units 
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Fig.  2.13  V(z)  of  Teflon  for  a  circular  aperture  lens:  (a)  observed  (solid),  (b)  calculated 
using  the  diffraction-corrected  ray  approach  (dashed)  and 
(c)  calculated  using  the  geometric  ray  theory  (dotted). 
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Fig.  2.14  Calculated  leaky  contribution  Vl(z)  for 


2.3.6  Calculation  of  the  total  transducer  response  V(z) 

To  calculate  the  final  transducer  response  V(z),  diffraction-corrected  components 
VG  and  Vl  are  added  in  accordance  with  (2.2.1).  The  resulting  V(z)  is  plotted  along 
with  the  experimentally  obtained  V(z)  in  Fig.  2.15,  and  with  the  Y(z)  calculated  using 
the  angular  spectrum  approach  in  Fig.  2.16.  Including  diffraction  effects  significantly 
improves  the  agreement  with  the  measured  result,  as  can  be  seen  by  comparing  Fig.  2.15 
with  Fig.  2.6,  which  displays  the  V(z)  calculated  while  ignoring  diffraction.  Further 
improvements  in  modelling  of  V(z)  near  z  =  0  would  require  the  corrections  in  Vl(z) 
mentioned  above. 
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using  the  diffraction-corrected  ray  approach  (dotted). 


Relative 

units 
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Fig.  2.16  V(z)  of  glass  for  a  circular  aperture  lens  calculated  using 

(a)  the  angular  spectrum  approach  (solid)  and 

(b)  the  diffraction-corrected  ray  approach  (dashed). 


2.4  Extraction  of  the  Rayleigh  wave  velocity  vr  from  the  V(z) 


It  has  been  shown  by  Bertoni  [2]  that  the  spacing  delta  z  between  the  V(z) 
minima  is  directly  related  to  the  Rayleigh  velocity  vr: 


a  -  ^*w 

Z  ”  2(l-COS  Or  ) 


(2.4.1) 


This  result  follows  from  the  ray  model  of  the  V(z)  [2]  and  is  valid  for  large  Izl  when  the 
detected  geometric  fields  have  nearly  planar  phase  fronts.  It  is  possible  to  deduce  vr 
from  the  V(z)  by  simply  measuring  the  separation  between  the  neighboring  minima  of 
V(z).  This  approach  relies  on  only  a  few  points  and  is,  therefore,  sensitive  to  the 
presence  of  noise  in  the  V(z).  An  alternative  method,  proposed  by  Kushibiki,  et.  al.  [2], 
finds  the  A  z  by  Fourier  transforming  the  pre-  processed  form  of  the  V(z)  to  obtain  the 
location  of  the  spatial  frequency  peak  corresponding  to  the  periodicity  in  the  ripple  in  the 
V(z).  This  method  has  the  advantage  of  using  all  of  the  collected  V(z)  data. 

The  following  is  a  partial  outline  of  an  improved  method  of  V(z)  processing, 
proposed  by  Briggs  [13],  which  will  be  implemented  here. 

The  microscope  used  in  this  work  measures  only  the  amplitude  and  not  the  phase 
of  the  transducer  voltage  V.  Thus,  when  V(z)  data  is  collected,  the  measured  quantity  is 

IV  (z)l  =  IVG(z)  +  VL(z)l 


The  square  of  the  amplitude  becomes 

IV(z)l2  =  IVGI2  +  2  iVGIIVLl  cos(argVG  -  argVL)  +  IVL!2  (2.4.2) 
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The  periodic  component  in  the  V(z)  is  due  to  the  cos(argVG  -  arg^)  factor  in  the  second 
term  which  is  the  result  of  the  interference  between  the  geometric  and  leaky 
contributions.  It  would  therefore  be  useful  to  isolate  this  factor  before  attempting  to 
determine  the  period  of  the  ripple  in  V(z)  since  the  other  terms  and  factors  would  obscure 
the  desired  peak  in  the  Fourier  spectrum. 

To  remove  IVGI  from  (2.4.2)  it  should  be  noted  that  the  functional  form  of  VG 
depends  only  on  the  lens  geometry  and  is  the  same  for  all  materials.  It  is  therefore 
possible  to  determine  Vq  by  measuring  the  V(z)  of  a  material  such  as  teflon,  which  has 
no  leaky  contribution.  The  measured  IV(z)l  of  teflon  must  then  be  properly  scaled  to 
account  for  the  difference  in  reflection  coefficients  of  glass  and  teflon.  One  practical 
way  to  do  so  with  measured  data  is  to  scale  the  IV(z)Ps  of  the  two  materials  to  match  at  a 
z  location  where  the  leaky  component  Vl(z)  is  negligible  or  zero,  such  as  for  z  «  0  or  z 
>0. 

Once  the  desired  IVq(z)I  of  glass  is  obtained,  a  reduced  form  of  IV(z)l,  called 
Vj(z),  can  be  found  from  (2.4.2),  where  Vj(z)  is  defined  as 

IV(z)l2  -  IVGI2  IVLI2 

VI(Z)  = - iy^j-^  +  2 1 VLI  cos(argVG  -  argVL)  +  .  (2.4.3) 

It  will  be  assumed  in  this  work  that  IVJ  is  a  slowly  varying  function  of  z  and  that  it  is 
small  compared  to  IVGI  so  that  there  is  no  significant  spectral  overlap  between  the 
cos(argVG  -  argVL)  term  and  the  other  terms  in  (2.4.3).  As  an  illustration,  a  Vj(z)  of 
glass  is  plotted  in  Fig.  2.17.  The  amplitude  of  the  Fourier  spectrum  of  Vi(z),  calculated 
by  a  Discrete  Fourier  Transform  (DFT)  routine,  is  displayed  in  Fig.  2.18.  The  location 
of  the  spectral  peak  in  Fig.  2.18  is  equal  to  1/Az.  The  Rayleigh  velocity  vR  for  glass, 
deduced  with  the  aid  of  (2.4.1),  is  found  to  be  3164  m/s.  When  either  the  angular 
spectrum  approach  or  the  ray  methods  are  used  to  theoretically  predict  the  V(z)  of  glass, 
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Fig.  2.17  Vj(z)  of  glass  obtained  from  measured  data  by 
V[(z)  =  ( V2(z)  -  V02(z)  )  /  VG(z). 


units 
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Fig.  2.18  Amplitude  of  Fourier  spectrum  of  Vj(z)  of  glass  calculated  by  DFT. 


the  best  fit  is  obtained  for  vR  =  3135  m/s.  which  is  assumed  to  be  the  nominal  value  of 
vR  for  this  particular  glass  sample.  Thus,  a  vR  has  been  deduced  from  an  experimental 
V(z)  with  an  error  of  about  1%. 

Applying  the  above  procedure  to  a  theoretically  generated  V(z)  gives  a  value  of 
vR  higher  by  about  1%  than  the  value  initially  assumed  during  the  generation,  which 
demonstrates  a  bias  consistent  with  the  processing  of  the  measured  V(z). 

Several  sources  of  this  error  have  been  found.  First  of  all,  as  mentioned  above, 
the  expression  (2.4.1]  for  the  spacing  between  minima  is  strictly  valid  only  for  z  «  0. 
For  z  near  0,  the  phase  of  the  Vq  exhibits  other  variations  with  z  which  alter  the  period 
in  the  interference  pattern  in  V(z).  Thus,  only  the  part  of  the  V(z)  far  from  the  focus  (z  = 
0)  should  be  used  in  the  extraction  of  vR.  Another  source  of  error  in  vr  is  the  residual  z- 
dependent  terms  in  (2.4.3),  like  IVl(z)I  and  IVl(z)I2/  IVq(z)I,  which  alter  the  spectrum  of 
Vj  and  should  therefore  be  removed  prior  to  Fourier  transformation.  Implementing  the 
above  corrections  in  the  processing  of  the  calculated  V(z)  returns  a  Rayleigh  velocity  that 
is  exactly  equal  to  the  initially  assumed  value. 

It  is  difficult  to  implement  these  corrections  for  the  processing  of  measured  data 
since  IVJ  is  not  a  measured  quantity  and,  therefore,  cannot  be  easily  removed  from  Vj. 
Furthermore,  as  IV(z)l  decreases  for  large  Izl,  the  signal- to-noise  ratio  deteriorates 
significantly,  making  the  V(z)  data  far  away  from  focus  unusable  for  transformation  into 
the  frequency  domain. 

In  summary,  the  above  method  can  be  used  to  find  the  Rayleigh  velocity  to  within 

1%. 
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Chapter  m 


Theoretical  Modelling  of  V(z)  for  the  Slot  Lens  Acoustic  Microscope 

Both  the  angular  spectrum  and  ray  approaches  are  used  to  model  the  slot  lens  for 
determining  V(z)  and  the  results  are  compared  with  experiment  for  several  materials. 
The  slot  lens  is  also  used  to  obtain  V(z)  of  Y-cut  quartz  for  a  number  of  slot  orientations 
relative  to  the  crystal  X-axis.  The  Rayleigh  wave  velocity  as  a  function  of  propagation 
direction  is  then  deduced  from  these  measurements  and  compared  with  accepted  values. 

3.1  Angular  Spectrum  Approach 

The  angular  spectrum  approach,  described  in  Chapter  II,  can  be  easily  adapted  for 
any  lens,  whose  aperture  plane  can  be  described  by  a  pupil  function  P(x,y).  Specifically, 
the  pupil  function  for  the  slot  lens  is  specified  by 

P(x,y)  =  circ  (^)  circ(^XR*y  - )  (3.1.1) 

where  w  is  the  half-width  of  the  slot  (w  <  Ra). 

To  find  V(z),  the  above  P(x,y)  is  substituted  into  expression  (2.1.2).  The 
resulting  V(z),  calculated  for  the  case  of  a  Teflon  sample  and  a  slot  width  w  =  1.5  mm  is 
plotted  in  Fig.  3.1  together  with  the  experimental  result.  A  V(z)  of  a  glass  sample  is 
compared  with  measurements  of  a  microscope  slide  in  Fig.  3.2. 

It  can  be  seen  that  calculated  and  measured  results  are  in  reasonable  agreement. 
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angular  spectrum  approach  (dashed). 


3.2  Ray  approach 


3.2.1  Calculation  of  the  specularly  reflected  contribution  Vq(z) 

The  slot  lens  pupil  function  P(x,y)  (3.1.1)  can  be  expressed  in  terms  of  cylindrical 
coordinates  (p,<J))  as 


P(p,<I>)  =  circ  circ  (j^) 


The  geometric  contribution  for  the  slot  lens  may  now  be  expressed  in  a  manner  similar  to 
expression  (2.3.14)  as 

Tlw(0)  Twl(O)  (R(0)  exp  (2i  kw  (D/n+f+z))  4 


VG(z)=- 


Mx 


7tRa2 


jt/2  Ra 


■J  /  exp(2S(^o)^iA,iexp(5l2)  pl(M^)P2(P)P<>P<it. 

0  0  * 

(3.2.1) 

where  Al,-  and  Bl,-  are  defined  in  (2.3.15).  In  (3.2.1)  integration  over  <j)  need  be  carried 
out  only  from  0  to  k/2,  owing  to  the  four-fold  symmetry  of  the  slot  aperture  lens.  The 
factor  of  4  is  then  introduced  to  account  for  the  contribution  of  all  four  quadrants. 

Using  the  slot  lens  pupil  function,  the  integral  over  p  in  (3.2.1)  can  be  carried  out 
analytically,  resulting  in: 


Tlw(0)  Twl(O)  m  exP  (2i  kw  (D/n+f+z)) 
VG(z)  = - - •  4 


k/2 


/  2|A1,t 


0 


nBlfl-u) 
ikwBlj2  -  2n(-u) 


*  [exp( 


&w  Pt2\  /~Pt2\ 

)  •  - 


2n(-u) 


1] 


d(J) . - 


(3.2.2) 
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where 


_  Wl 
”  COS  (j)  ’ 

<|>  <  COS-1  (^) 

EL 

pT  =  Ra  > 

<t>  >  cos-1  (j^)  . 

and 

wl  =wMj, 

MX<1 

ii 

>  1. 

Introduction  of  the  slot  aperture  reduces  the  symmetry  of  the  lens,  thus  the  above 
integral  cannot  be  further  evaluated  analytically,  and  numerical  integration  is  used 
instead.  The  result  of  this  calculation  is  plotted  in  Fig.  3.3  and  compared  with 
experimental  measurements.  The  ray-theoretical  Vq(z)  does  not  agree  with 
measurement  as  well  as  the  angular  spectrum  Vq(z)  (Fig.  3.1).  This  could  be  due,  in 
part,  to  the  fact  that  the  geometric  approach  ignores  diffraction  in  water,  which  may  be 
more  significant  for  the  slot  aperture  lens  than  for  circular  lens. 

The  advantage  of  the  ray  approach  outlined  above  is  that  it  requires  only  a  single 
integration  over  <|>.  The  angular  spectrum  approach  requires  a  double  integral  over  (x,y), 
which  is  much  more  computationally  intensive. 

3.2.2  Calculation  of  the  leaky  contribution  Vl(z) 

As  discussed  in  Chapter  II,  leaky  rays  emanating  from  the  sample  surface  form  a 
ring  focus  of  radius  pR  in  the  lens  rod  of  the  circular  lens.  The  slot  aperture  allows  only 
portions  of  this  ring  to  be  formed.  Thus,  neglecting  diffraction,  the  leaky  contribution 
Vl(z)  in  the  slot  lens  is  the  same  as  for  circular  lens  but  scaled  by  the  ratio  of  the  length 
of  arc  exposed  by  the  slot  to  the  total  circumference  of  the  focal  ring  or 
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VL(z)  =  Vhtirc  (z)  * 


(3.2.3) 


2  sin-1  ( w  /  Qp  ) 
TC 


where  w  is  the  half-width  of  the  slot. 

3.2.3  The  total  voltage  V(z) 

The  total  voltage  response  V(z)  of  the  slot  aperture  lens  for  a  glass  sample  is 
obtained  by  adding  Vq(z)  and  Vl(z)  and  plotted  in  Fig.  3.4  for  comparison  with  the 
measurement. 
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the  diffraction-corrected  ray  approach  (dashed). 
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A 


diffraction-corrected  ray  approach  (dashed). 


3 3  Extracting  vR(<|))  of  an  anisotropic  material  using  the  slot  lens 

A  slot  lens  with  a  slot  width  w  of  1  mm  was  constructed  and  used  to  obtain  V(z) 
measurements  of  a  Y-cut  quartz  sample  for  various  orientation  of  the  slot  relative  to  the 
crystalline  X-axis.  The  resulting  V(z)  curves  for  orientation  angles  of  0°  and  90°  are 
displayed  in  Figs.  3.5  and  3.6,  respectively.  To  obtain  the  Rayleigh  wave  velocity  vR  as 
a  function  of  propagation  direction  <j),  these  curves  were  processed  using  the  algorithm 
described  in  Chapter  n.  The  resulting  values  are  plotted  in  Fig.  3.7  along  with  the 
accepted  dispersion  curve  vR(<|))  (Kushibiki  [4]).  For  illustration,  a  processed  Vj(z) 
curve,  defined  in  Chapter  II,  is  shown  in  Fig.  3.8  for  <j)  =  0°.  Higher  precision  in 
calculating  vR  from  experimental  curves  require  more  periods  in  Vi(z).  This  occurs  in 
materials  with  low  Rayleigh  wave  velocities,  since  the  period  of  the  interference 
decreases  with  decreasing  vR,  as  can  be  seen  from  (2.4.1).  For  materials  with  high 
Rayleigh  velocities  higher  acoustic  wave  frequencies  should  be  used.  Alternatively,  it  is 
possible  to  simply  measure  V(z)  over  a  greater  range  in  z,  but  this  would  require  a  higher 
signal-to-noise  ratio  and  possibly  a  different  lens  geometry  since  the  signal  level  drops 
off  significantly  for  large  defocus  distance  Izl. 

It  has  been  shown  that  a  slot  lens  does  indeed  launch  Rayleigh  waves  in  a 
preferential  direction  and  can  be  used  to  obtain  Rayleigh  wave  velocities  as  a  function  of 
propagation  direction  in  anisotropic  crystals. 
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Fig.  3.5  Observed  V(z)  of  Y-cut  quartz  with  slot  oriented  along  the  crystalline  x-axis 
(♦  =  0). 
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Fig.  3.6  Observed  V(z)  of  Y-cut  quartz  with  slot  oriented  at  90°  to  the  crystalline  x-axis. 


<(>  (degrees) 


Fig.  3.7  Rayleigh  wave  velocity  vR  of  Y-cut  quartz  as  a  function  of  direction  of 

propagation  <j>  measured  from  the  crystalline  x-axis:  accepted  (solid  line)  and 
measured  using  the  slot  lens  (open  boxes). 
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Fig.  3.8  Vj(z)  of  Y-cut  quartz  obtained  from  data  measured  with  slot  oriented  along  the 
crystalline  x-axis  (<J>  =  0). 


Chapter  IV 

Modelling  of  V(x)  measurements  across  a  surface  breaking  crack 

using  ray  theory. 

This  chapter  is  concerned  with  a  theoretical  treatment  of  the  microscope  response 
of  a  sample  surface  containing  a  discontinuity  such  as  a  surface-breaking  crack.  This 
defect  will  be  assumed  to  be  of  infinititesimal  thickness,  straight,  symmetrical  in  a  planar 
surface  and  long  compared  with  the  diameter  of  the  spot  produced  by  the  incident 
acoustic  beam.  Of  particular  interest  are  line-scan  measurements,  called  V(x),  obtained 
by  keeping  the  lens  at  a  constant  defocus  distance  z  and  scanning  in  x  perpendicular  to 
the  discontinuity  (see  Fig.  1.5).  Crack  displacement  x  is  defined  as  the  distance  between 
the  lens  axis  and  the  crack,  with  x  =  0  corresponding  to  the  lens  being  centered  directly 
above  the  crack.  The  analysis  will  first  involve  the  slot  lens,  where  the  angular  range  of 
surface  ray  paths  is  restricted,  and  then  the  circular  lens  can  be  viewed  as  simply  a 
special  case  of  the  slot  lens  with  slot  half-width  w  equal  to  the  aperture  radius  w  =  Ra. 

The  total  response  V(x,z)  of  the  microscope  to  a  surface  containing  a  crack  may 
be  expressed  as 


V(x,z)  =  VG(z)  +  VT(x,z)  +  Vrefl(x,z),  (4.1) 

where  Vq(z)  is  the  specularly  reflected  component,  Vj(x,z)  is  due  to  the  leaky  waves 
transmitted  through  the  crack  as  well  as  those  leaky  waves  that  propagate  unimpeded  by 
the  crack,  and  Vrefl  (x,z)  is  the  contribution  of  the  leaky  waves  scattered  from  the  crack. 
The  ray  paths  involved  in  determining  these  contributions  are  shown  in  Fig.  4.1.  Here 
both  the  crack  and  the  lens  are  assumed  to  have  inversion  symmetry,  V(x,z)  =  V(  -x,z), 
i.e.  the  same  response  is  detected  on  either  side  of  the  crack.  It  will  be  assumed  that  when 
a  surface  wave  strikes  a  crack,  a  portion  of  the  incident  wave  is  transmitted  with  a  field 
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Fig.  4.1  Surface-wave  and  specular  ray  paths  for  various  contributions  to  V(x,z): 
(a)  transmitted  through  the  crack  ,  (b)  reflected  from  the  crack,  and 
(c)  specularly  reflected  from  the  surface. 


transmission  coefficient  Tc,  and  a  portion  is  reflected  with  a  field  reflection  coefficient 
Rc.  The  dependence  of  Tc  and  R<,  on  the  angle  of  incidence  of  the  surface  wave  will  be 
ignored  on  the  basis  that  the  predicted  variation  near  normal  incidence  is  weakly 
dependent  on  angle,  as  shown  by  Angel  and  Achenbach  in  [3]. 


4.1  Specularly  reflected  component  Vq(z) 

The  crack  is  assumed  to  be  of  infinitesimal  thickness,  therefore  the  specularly 
reflected  component  Vq  is  unperturbed  by  the  crack.  Since  Vq  is  the  same  as  for  the 
defect-free  surface,  it  is  independent  of  the  crack  position  x  and  only  depends  on  defocus 
distance  z.  Thus,  Vq(z)  is  evaluated  by  the  methods  described  in  Chapters  II  and  III  that 
were  developed  for  defect-free  substrates. 

4.2  Transmitted  component  Vj(x,z) 

It  is  instructive  first  to  review  the  leaky  ray  structure  for  the  defect-free  surface 
when  the  lens  is  closer  than  focus  or  z  <  0.  Some  rays  cross  the  lens  surface  at  p  =  pR 
and  strike  the  sample  surface  at  the  Rayleigh  critical  angle  0R,  thus  launching  Rayleigh 
surface  waves  in  the  plane  of  incidence.  The  locus  of  all  such  launch  points  forms  a 
circle,  called  Rayleigh  circle,  in  the  case  of  a  circular  lens,  or  restricted  circular  arcs  in 
the  case  of  the  slot  lens.  The  radius  c  of  the  Rayleigh  circle  is  found  from 

c  =  Izl  tan  (0R) .  (4.2.1) 

The  surface  waves  focus  at  the  intersection  with  the  lens  axis  and  diverge  beyond.  As 
discussed  in  Chapter  II,  the  leaky  rays  form  a  circle-shaped  focus  of  radius  pR  within  the 
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lens  rod,  located  at  fx  above  the  lens  surface,  which  is  limited  to  circular  arcs  for  the  slot 
lens.  The  leaky-ray  bundle  reaching  the  transducer  is  defined  by  rays  that  cross  the  lens 
surface  from  p  =  Pr^  to  p  =  pRmax*  as  shown  in  Fig.  4.2.  This  ray  bundle  originates 
from  points  on  the  sample  surface  in  the  radial  plane  lying  between  p  =  cmin  and  p  = 
Cjnax-  The  geometric  quantities  cmin,  cmax,  pmin,  pmax  are  derived  in  Appendix  A. 
Thus,  the  contribution  Vl(z)  is  due  to  a  distribution  of  effective  leaky  ray  sources  on  the 
sample  surface,  whose  locus,  shown  in  Fig.  4.3,  is  two  sections  of  a  ring,  bounded  by 
circles  of  radii  cmin  and  cmax  and  the  projection  of  the  slot  aperture  onto  the  sample 
surface. 

When  the  crack  is  sufficiently  displaced  from  the  lens  axis,  some  surface  rays  are 
leaked  into  water  without  interacting  with  the  crack,  and  therefore  their  contribution  to 
the  transducer  voltage  is  the  same  as  in  the  defect-free  case.  Voltage  due  to  those  rays 
that  are  transmitted  through  the  crack  is  scaled  by  the  crack  transmission  coefficient  Tc 
(see  Fig.  4.4). 

Transmitted  contribution  Vj(x,z)  is  given  by 

VT  =  VL(z)  •  (Qr(x,z)  +  Ql(x,z)),  (4.2.2) 

where  VL(z)  is  the  contribution  to  the  V(z)  due  to  the  leaky  surface  waves  for  the  defect- 
free  surface,  and  Qr(x,z)  and  Ql(x,z)  are  fractional  contributions  to  V(x,z)  from  the  rays 
re-radiated  into  water  on  the  right  and  left  sides  of  the  lens,  respectively,  as  shown  in  Fig. 
4.1a.  In  the  absence  of  a  crack  Qr  =  Ql  =  0.5. 

First,  only  rays  launched  on  the  left  side  of  the  crack  and  propagating  to  the  right 
will  be  considered.  The  VT(x,z)  contribution  to  the  V(x,z)  will  be  assumed  to  be 
proportional  to  the  total  field  flux  leaked  into  water  from  the  area  shown  shaded  in  Fig. 
4.3.  This  assumption  is  most  valid  for  a  long  lens  rod  where  the  phase  variation  of  the 
field  at  the  transducer  is  small.  The  factor  Qr  can  be  seen  to  be 
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Fig.  4.2  Leaky  ray  bundle  in  the  x-z  plane. 
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Fig.  4.3  Distribution  of  the  leaky  ray  flux  (shaded)  that  excites  the  transducer,  defined 
on  the  sample  surface.  Note  also  the  definition  of  <(>i  and  <])max  in  the  sample 
plane. 
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Fig.  4.4  Geometry  of  a  surface  ray  is  excited  on  the  left  and  leaks  energy  to  the  right  of 
lens  axis  .  After  the  ray  crosses  the  crack  the  leaky  field  is  reduced  by  the 
factor  Tc. 


_  QiT  +  1-Qi 

Qr  -  2 


(4.2.3) 


where  Qj  is  the  fraction  of  the  ray  flux  intercepted  by  the  crack  and  therefore  scaled  by 
Tc,  and  (1  -  Ch)  is  the  fraction  of  the  ray  flux  that  is  unobstructed.  It  can  be  seen  from 
Fig.  4.4  that  the  factor  Qj  can  be  expressed  as: 


Ql 


Area  of  obstructed  part  of  Rayleigh  ring 
Total  area  of  Rayleigh  ring  section 


(4.2.4) 


The  above  relation  simply  views  all  leaky  launch  points  as  equally  contributing  to 
Vj(x,z).  Total  area  of  the  Rayleigh  ring  may  be  found  approximately  from  the  relation 


24*1  ft(cmax^  ~  cmin^) 
n  2 


~  4>1  (cmax^  "  cmin^)- 


[4.2.5) 


where  (jjj,  illustrated  in  Fig.  4.3  is  the  effective  half-angle  of  the  ring  arc,  defined  in 
terms  of  the  slot  half-width  w  at  the  lens  aperture  by 


<>!  =  sin*1  (w  /  pR), 

=  K/2, 


W<pR 

W>pR. 


(4.2.6) 


The  second  equality  in  (4.2.6)  applies  to  situation  where  none  of  the  area  defined  by  cmin 
<  pR  <  cmax  is  obstructed  by  the  slot.  All  azimuthal  angles  <J)  are  defined  relative  to  the 
direction  of  the  x-axis.  As  the  lens  is  translated  to  the  left  along  the  x-axis,  there  are  four 
distinct  ranges  of  crack  position  x,  shown  in  Fig.  4.5.  In  Fig.  4.5  <|>max  is  the  azimuthal 
angle  coordinate  which  forms  one  of  the  boundaries  of  the  Rayleigh  ring  arc,  whose 
value  when  defined  at  the  lens  aperture  is 
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(C) 


(d) 


Fig.  4.5  Distinct  ranges  of  crack  displacement  that  lead  to  different  geometric  relations 
in  determining  Qr: 

(a)  0  <  X  <  Cmjn  ,  (b)  Cjjjjjj  cos  <|>max  <x<  cmin 

(c)  cmtn  <  x  <  cmax ,  (d)  x  >  cmax . 
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w  <  Pmin 
w  >  Pmin 


(4.2.7) 


4>max  =  sin'1  (w  /  pmin). 


=  k/2. 


For  Range  1,  where  0  <  x  <  cmin  cos  <J>max ,  all  surface  rays  are  intercepted  by  the 
crack  before  leaking  into  water  and,  therefore,  Qj  =  1.  Substituting  Ql  into  (4.2.3)  gives 
Qr  =  Tc/2. 

For  Range  2,  where  cmin  cos  <J>max  <  x  <  cmin,  some  ray  paths  are  no  longer 
obstructed  by  the  crack.  Upon  examination  of  Fig.  4.5b  and  in  accordance  with  (4.2.4), 
Ql  may  be  found  approximately  from  the  expression 


cmax  cos4(x/p) 

Ql  =  — J— — 2\  J  J  d(j)  p  dp  = 

4>1  (cmaxZ  -  cmin  ) 


cmin  0 


^max 


01  (cmax^  *  cmin^) 


J  p COS'^x/p)  dp  . 


'irnn 


(4.2.8) 


Similarly,  for  Range  3,  where  cmjn  <  x  <  cmax,  it  follows  from  Fig.  4.5c  and 
(4.2.4)  that 


-max 


Ql  = 


4>1  (c, 


max 


^min 


•  J  p  cos’^x/p)  dp  . 


[4.2.9) 


For  Range  4,  where  x  >  cmax,  all  surface  rays  propagate  unobstructed  by  the 
crack,  thus  Qi  =  0,  and  Qr  =  0.5. 
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Surface  rays  that  are  excited  on  the  right  and  detected  on  the  left  will  be  now 
considered.  As  can  be  seen  in  Fig.  4.6  there  are  three  distinct  ranges  of  crack  position  x. 

For  Range  1,  where  x  <  c  cos  <(»|_,  all  rays  are  intercepted  by  the  crack,  so  that  Ql 
=  Tc/2. 

For  Range  2,  where  c  cos  <  x  <  c,  it  can  be  seen  from  Fig.  4.6b  that 


Ql  =  0.5 


2  4>T  +  2($H>) 
2<h 


<!>  T  +  <|>i  -  <)> 
2<h  ’ 


(4.2.10) 


where  <{>  in  Fig.  4.6b  is  defined  as  <j>  =  cos-1  (x  /  c). 

For  Range  3,  where  x  >  c,  all  rays  pass  unobstructed,  so  that  Ql  =  0.5. 

The  transmitted  component  Vj(x,z)  can  be  viewed  as  a  reduction  of  the  leaky 
component  Vl(z)  due  to  the  presence  of  the  crack.  The  above  analysis  enables  the 
determination  of  Vt(x,z)  for  z  <  0.  When  z  >  0,  the  surface  waves  diverge  away  from 
the  lens  and,  thus,  there  is  no  leaky  component  Vl(z)  intercepted  by  the  transducer,  as 
discussed  in  Chapter  n.  Consequently,  Vt(x,z)  =  0  for  z  >  0. 
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(  C  ) 


Fig.  4.6  Distinct  ranges  of  crack  displacement  that  lead  to  different  geometric  relations 
in  determining  Ql:  (a)0<x  <  c-cos  ,  (b)  c-cos  (J)j  <  x  <c  ,  (c)  x>c 
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4.3  Reflected  surface  wave  contribution  VreflW: 


In  Fig.  4.7,  the  reflected  contribution  to  V(x,z)  has  two  components,  defined  by 

Vrefl(x,z)  =  Vrr(x,z)  +  Vrl(x,z)  .  (4.3.1) 

where  Vrr  is  the  contribution  due  to  the  surface  waves  incident  on  the  discontinuity 
from  the  right  (Fig.  4.7a)  and  Vrl  is  due  to  the  waves  incident  from  the  left  (Fig.  4.7b). 
Further  examination  of  the  figure  reveals  that  Vrr  =  0  for  x  >  c  and  z  <  0.  In  this  case 
the  only  contribution  to  Vrefl  is  from  the  left.  Let  Vrefl(x,z)  =  VR(x,z)  in  this  range. 
First,  Vr(x,z)  will  be  found  and  then  both  VRR(x,z)  and  VRL(x,z)  will  be  expressed  in 
terms  of  VR(x,z)  for  all  x  and  z. 

The  various  field  amplitude  factors  due  to  focusing  and  phase  delays  due  to 
propagation  are  now  to  be  determined.  Recall  that  in  Chapter  II  determination  of  the 
leaky  component  Vl(z),  a  ray  field  of  amplitude  IU(pR)l  is  incident  on  the  lens  surface  at 
the  critical  radius  p  =  pR  where  it  is  refracted  into  the  water  to  strike  the  sample  surface 
at  p  =  c.  Spherical  focusing  by  the  lens  and  propagation  in  the  lens  rod  and  water 
introduces  the  amplitude  factor  p  /  c  and  the  phase  factor  exp  i  kw  (D/n  +  f  +  z  sec  0R)  ( 
see  Fig.  4.2).  Propagation  along  the  surface  to  the  crack  results  in  the  phase  factor  exp 
ikp(x  -  z  tan  0R),  where  kp  is  the  pole  of  the  water-substrate  reflection  coefficient:  kp  = 
kR  +  ictL  ([2]).  Reflection  from  the  crack  introduces  the  factor  Rc.  Traversing  the  return 
path  to  the  transducer  causes  the  repetition  of  the  phase  changes  involved  in  forward 
propagation  as  well  as  an  additional  amplitude  term  exp  -ocl  (x  -  z  tan0R)  due  to  leakage 
into  water.  Transmission  through  the  lens-water  interface  is  taken  into  account  by  using 
the  field  transmission  coefficients  Tlw(Pr)  and  T^l(pr)  mentioned  in  Chapter  II.  In 
addition,  as  shown  by  Bertoni  [2],  there  is  a  coefficient  E  associated  with  the  excitation 
and  re-radiation  of  the  surface  waves,  given  by 
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(  a) 


(b) 


Fig.  4.7  Focusing  of  surface  rays  reflected  from  the  crack  on  the  (a)  left  and 
(b)  right  side  of  the  Rayleigh  circle. 
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(4.3.2) 


/  2  71  Izl  ,  . 

E='2“L  V 

When  the  distance  between  the  excitation  and  re-radiation  of  surface  waves  is  small:  (x- 
Xj)  <  'sjlzl  7^  /  cos3  0R,  this  coefficient  must  be  modified,  as  described  by  Bertoni  [2]. 
Usually,  in  modelling  surface  wave-crack  interaction,  such  correction  will  become 
necessary  only  when  the  crack  is  near  the  edge  of  the  Rayleigh  circle,  i.e.  x  =  c.  In  this 
region,  the  Rayleigh  wave  is  launched,  scattered  from  the  crack,  and  re-radiated  within  a 
very  short  distance.  In  this  work  expression  (4.3.2)  will  be  used  to  determine  E  for  all  x. 

The  focusing  and  spreading  of  surface  wave  rays  for  several  Rayleigh  wave  ray 
paths  is  depicted  in  Fig.  4.8.  The  surface  waves  pass  through  a  focus  where  the  lens  axis 
intersects  the  sample  surface  and  then  strike  the  crack,  located  at  a  distance  x  from  the 
lens  axis,  which  acts  as  a  plane  mirror  with  a  field  reflection  coefficient  Rc.  The  surface 
waves  reflected  from  the  crack  appear  to  come  from  a  point  source  located  at  a  distance  2 
x  from  the  lens  axis  denoted  by  the  point  F  in  the  Fig.  4.8.  This  point  source  is  thus 
located  at  the  mirror  image  of  the  point  defined  by  the  lens  axis  intersection.  These 
reflected  waves  propagate  to  the  left  until  they  reach  the  neighborhood  of  the  Rayleigh 
circle,  where  the  ray  fields  that  leak  into  the  water  are  refracted  by  the  lens  and  strike  the 
transducer.  Since  the  surface  waves  are  launched  with  a  radius  of  curvature  of  c  and  are 
re-  radiated  with  a  radius  of  curvature  c  +  2x,  the  field  amplitude  is  modified  by  a  factor 
sjc  1  (c+2x).  In  addition,  the  planar  focus  at  the  lens  axis  introduces  a  phase  change  of 
-n/2. 

Each  of  the  surface  wave  rays  (  e.g.  FB,  FA,  and  FC  in  Fig.  4.8)  radiate  into 
water  thereby  creating  a  parallel  family  of  rays  at  angle  0R  with  respect  to  the  normal  to 
the  surface.  Since  parallel  rays  do  not  spread,  there  is  only  angular  spreading  along  the 
arc  BAC.  When  the  radiated  ray  bundle  reaches  the  lens,  it  will  have  a  radius  of 
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curvature  of  pR+2x,  and  thus  a  change  of  amplitude  of  -\y(c+2x)/(pR+2x)  resulting  from 
ray  spreading  in  the  water. 

Since  the  bundle  of  radiated  rays  spreads  only  in  the  angular  direction,  passage 
through  the  spherical  lens  results  in  two  astigmatic  foci.  The  rays  radiated  along  FB,  FA, 
and  FC  in  Fig.  4.8  are  parallel  to  each  other  and  therefore  focus  at  fx  above  the  lens 
surface.  For  the  direction  FA  the  focal  length  fx  is  the  same  as  fx  described  in  Chapter  II 
(see  Fig.  2.5).  For  the  focusing  of  leaky  rays  in  the  case  of  the  defect- free  surface 


(4.3.3) 


See  Appendix  A  for  derivation  of  (4.3.3).  It  is  assumed  here  that  fx  is  the  same  for  all 
directions.  This  approximation  is  most  valid  for  ray  paths  close  to  FA,  i.e.  nearly  normal 
to  the  crack.  The  rays  that  strike  the  crack  at  large  oblique  angles  are  strongly  refracted 
by  the  lens  and  are  likely  to  miss  the  transducer.  For  the  slot  lens  all  rays  incident  on  the 
crack  strike  at  nearly  normal  incidence. 

Those  rays  that  were  are  radiated  along  arcs  parallel  to  BAC  in  Fig.  4.8  appear  to 
originate  from  a  point  on  a  vertical  line  going  through  F,  shown  in  Fig.  4.9.  The  lens 
forms  an  image  of  this  point  source  at  fy  above  the  lens  surface,  where  fy  is  shown  in 
Appendix  A  to  be 


For  the  lens  used  in  this  work  and  for  x  >  c  it  follows  that  0  <  fy  <  D  and,  thus,  there  is  a 
focus  in  the  lens  rod.  The  ray  field  therefore  suffers  a  phase  change  of  -it/2  due  to 
propagation  through  each  of  the  two  foci  at  fx  and  fy  from  the  lens  surface  as  well  as 
additional  factors  of  /  (D-fx)  and  -\Jfy  /  ID-fyl  due  to  the  spreading  in  the  lens  rod. 


78 


1 


1 


Fig.  4.9  Focus  formed  in  the  lens  rod  by  the  rays  re-radiated  along  the  arc  B  AC 
in  Fig.  4.8. 


79 


These  factors  may  be  approximated  using  fx«D  and  fx  —  f/n  by 


Vy(w,)“i/yo 

-y/fy/ID-fyl-^/fc/nlD-fyl) 

Finally,  the  astigmatic  beam  illumination  of  the  transducer  has  a  phase  curvature  of 

&w  (x~Pr)2  &wy2 

exp  2nD  *  exp2n(D-fy) 


Combining  and  simplifying  algebraically  all  of  the  above  factors  gives  the  field, 
illuminating  the  transducer  as 

Rc  TLw(Pr)  tWl(Pr)  E  exp  2i(kwD/n+kw 
•  exp  2i  kp(x  -  z  tan  0r)  •  exp  ( -i  3tc/2) 

PR  /  D  ftwfrl-PR)2 

'  D^/2^Jk'^/D-fy  'eXp  2nD 

where  the  function 

Gate(y)  =1,  lyl  <  1 

=  0 ,  otherwise. 

takes  into  account  the  spatial  bounds  imposed  on  the  beam  by  the  lens  aperture.  Using 
the  divergence  of  the  beam  to  project  the  aperture  onto  the  transducer  plane  (Fig.  4.10), 
ymax  is  found  to  be 


f+z  sec  0r) 

a/? 

fcwyi2 


exp 


Gate 


2n(D-fy)  —  Vymax 


(4.3.5) 
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Fig.  4.10  Illumination  of  the  transducer  by  a  bundle  of  leaky  rays,  shown  in  y-z  plane, 
reflected  from  the  crack. 


81 


Ymax  “ 


(4.3.6) 


ID-fyl 

a_i ijr 


The  aperture  radius  Ra  is  equal  to  transducer  radius  Rj  for  the  circular  lens  used  in  this 
study.  It  follows  from  (4.3.6)  that  ymax  >  Rj  and  the  transducer  is  fully  illuminated  for 
fy  <  D/2.  To  find  ymax  for  the  slot  lens,  the  half-width  w  of  the  slot  should  be  used  in 
(4.3.6)  instead  of  Ra.  It  then  follows  that  in  the  case  of  the  slot  lens,  the  transducer  is 
fully  illuminated  for  fy  <  D/(Rj/w  +  1). 

Integrating  the  field  (4.3.5)  over  the  transducer  and  normalizing  to  the  transducer 
area,  the  voltage  response  VR  is  found  to  be 


Vr(x)=Rc  TLW(pR)  TWL(pR)  E  exp  2i(kwD/n+kw  f+z  sec  0R) 

jf  fx  pR 

•  exp  2i  kp(x  -  z  tan  0R)  •  exp(  -i3ic/2) 


K!(x),  (4.3.7) 


where 

1  f  /  D  ikw(xrpR)2  ikw  yj2 

K‘(X)=  5*5  '  J  exP  2nD  ’  exP2^(DZ)dxld>'l  (43.8) 

transducer  v  3  . 

The  x-dependent  factor  in  (4.3.7)  is  Kj(x)  exp  (2ikpx)  /  \J\x\ . 

As  Kl(x)  is  not  a  function  of  z  it  can  be  pre-calculated  once  and  used  for  all  z  in 
obtaining  numerical  solutions.  The  amplitude  and  phase  of  Kl(x)  were  evaluated 
numerically  and  plotted  in  Figs.  4.11  and  4.12  for  the  case  of  a  glass  substrate  and 
dimensions  of  the  lens  used  in  the  measurements.  Since  both  amplitude  and  phase  of  Kj 
are  nearly  constant  except  for  small  x,  the  value  of  Kj(x)  at  any  x  sufficiently  far  from 
the  origin  can  be  evaluated  just  once  for  a  particular  lens  and  used  for  all  large  values  of 
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Fig.  4.11a  Amplitude  of  Kj(x)  calculated  for  a  circular  lens. 
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Fig.  4.lib  Phase  of  Kj(x)  calculated  for  a  circular  lens. 
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Fig.  4.12b  Phase  of  Ki(x)  calculated  for  a  slot  lens. 


x.  It  follows  that  V(x)  varies  as  exp(2ikpx)  /  -s/ixi  for  large  x. 
At  fy  =  D,  which  corresponds  to 

(f-R)pR 
x  ”  2(D-fx)  ’ 


(4.3.9) 


the  focal  point  lies  in  the  plane  of  the  transducer,  resulting  in  rapid  phase  variations  of 
the  integrand  in  (4.3.8).  At  this  singular  point  geometric  ray  approximations  fail  and 
diffraction  effects  must  be  considered.  However,  as  is  evident  from  (4.3.4),  fy  is  a 
rapidly  varying  function  of  x  (for  small  x)  and,  thus,  the  approximations  quickly  become 
reasonable  only  a  small  distance  from  this  critical  point. 


The  reflected  contribution  Vren(x,z)  for  z  <  0 

The  entire  reflected  voltage  contribution  Vrefl(x,z)  will  now  be  considered  when 
the  crack  lies  within  the  Rayleigh  circle,  i.e.  z  <  0  and  x  <  c,  shown  in  Fig.  4.7a  .  The 
rays  incident  from  the  left  give  exactly  the  same  expression  as  found  for  Vr(x)  except 
that  for  fy  >  D  there  is  only  one  astigmatic  focus  in  the  lens  rod,  which  corresponds  to 
0  <  x  <  2(D-fx)  /  (f-R)  pR.  In  this  case  the  ray  field  does  not  suffer  a  -7t/2  phase  shift  due 
to  passing  through  the  focus  in  the  corresponding  range  of  x,  or  explicitly 

Vrl(x)  =  VR(x)  exp  ( irc/2),  0  <  x  <  2(D-fx)  /  [(f-R)  pR] 

Vrl(x)  =  Vr(x),  otherwise.  (4.3.10) 


To  find  the  contribution  due  to  rays  incident  on  the  right,  -x  should  replace  x  in  the 
expression  (4.3.7).  As  shown  in  Fig.  4.7b,  the  Rayleigh  waves  are  reflected  from  the 
crack  before  reaching  the  focus  at  the  lens  axis.  The  reflected  rays  will  still  converge  to 
a  focus  located  a  distance  2x  from  the  lens  axis.  If  2x  >  c,  the  rays  will  leak  into  water 
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and  focus  above  the  sample  surface.  In  either  case,  a  real  focus  is  formed  and  the 
previous  expression  derived  for  VR(x,z)  applies.  It  should  also  be  noted  that  since  fy  <  0 
for  Ixl  small  and  x  <  0,  there  is  only  one  real  focus  in  the  lens  rod  and  therefore: 

Vrr(x,z)  =  VR(  -x,z)  exp(  (ire/2).  (4.3. 1 1) 

As  Ixl  increases,  fy  for  VrrCx.z)  determination  will  eventually  become  positive, 
however,  in  this  case  the  crack  will  be  outside  the  Rayleigh  circle  for  practical  defocus 
distances  z  and  moreover  VrrCx.z)  =  0  for  Ixl  >  c. 

The  reflected  contribution  Vren(x,z)  for  z  >  0 

It  can  be  seen  from  Fig.  4.13  that 

Vrj^x.z)  =  0 ,  for  z  >  0,  x  >  0 

Vrl/x.z)  =  0  for  z  >  0,  Ixl  <  c.  (4.3.12) 

For  z  >  0  and  x  >  c,  the  diverging  beam  incident  on  the  sample  surface  has  passed 
through  the  spherical  lens  focus  and  has  therefore  suffered  a  phase  shift  of  -n.  However, 
there  is  no  real  Rayleigh  wave  focus  on  the  sample  surface  in  this  case  since  the  incident 
beam  is  diverging  when  it  strikes  the  sample  surface  and  the  phase  shift  of  -re/2  included 
in  derivation  of  VR(x,z)  no  longer  applies.  Thus,  for  z  >  0 

Vrefl(x,z)  =  VL(x,z)  =  VR(x,z)  exp(-bt/2),  x  >  c 

=  0,  x<c.  (4.3.13) 
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Fig.  4.13  Surface  wave  interaction  with  a  crack  for  z  >  0 . 
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4.4  The  total  voltage  response  V(x,z) 


The  total  transducer  voltage  response  function  V(x,z)  is  found  by  adding  the 
specularly  reflected,  transmitted  and  crack-scattered  contributions  as  per  expression 
(4.1).  Typical  V(x,z),  calculated  for  circular  and  slot  lenses  for  defocus  distance  z  =  -400 
pm  and  crack-lens  displacement  0  <  x  <  500  |im,  are  plotted  in  Figs.  4.14  and  4.15  , 
respectively. 

As  mentioned  earlier,  it  is  assumed  here  that  the  surface  wave  reflection  and 
transmission  coefficients  of  the  crack,  R<.  and  Tc>  respectively,  are  constant  for  all  angles 
of  incidence.  It  is  therefore  necessary  to  use  proper  effective  values  for  Rc  and  Tc  . 
Angel  and  Achenbach  [3]  have  carried  out  the  analysis  to  find  the  scattering  coefficients 
for  various  angles  of  incidence  and  crack  depths.  Some  of  their  results  are  plotted  in  Fig. 
4.16.  The  crack  used  in  this  work  was  very  deep  and,  thus,  the  crack  depth  was  assumed 
to  be  infinite.  Since  the  circular  lens  excites  surface  waves  that  strike  the  crack  over  a 
large  range  of  incident  angles,  it  was  decided  to  use  =  1  and  Tc  =  0  as  the  effective 
values  in  the  calculations  for  this  lens.  This  assumption  is  most  valid  for  small  x.  As  x 
increases,  the  rays  incident  at  oblique  angles  that  have  a  large  miss  the  transducer  so 
only  rays  that  are  nearly  normal  with  lower  effective  R^  and  larger  effective  Tc  are 
detected. 

The  surface  waves  launched  by  the  slot  lens  propagate  in  a  narrower  range  of 
directions  where  the  largest  angle  of  incidence  on  the  crack  is  given  by  tan*1  (w/pR).  For 
the  lens  used  here,  the  above  angle  is  22°.  The  effective  values  for  the  slot  lens  are 
therefore  chosen  as  Rc  =  .54  and  Tc  =  0,  which  are  strictly  valid  only  for  normally 
incident  surface  waves. 
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Fig.  4.14  V(  x,  z  =  -400  |im)  calculated  for  a  circular  lens. 
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Chapter  V 

Experimental  measurement  of  V(x) 


In  this  chapter  a  comparison  will  be  made  between  experimentally  measured  V (x) 
and  V(x)  calculated  using  the  ray  theory  discussed  in  Chapter  IV.  An  algorithm  for 
extracting  the  reflection  coefficient  of  the  crack  from  the  observed  V(x)  data  will  be 
presented  and  used  to  determine  the  Rc  value  for  a  cracked  glass  sample. 

5.1  Experimental  measurement  of  V(x) 

To  obtain  a  suitable  sample,  1  mm  Coming  glass  microscope  slides  were  stressed 
with  a  pointed  object  to  introduce  a  linear  localized  defect.  The  sample  was  then  heated 
on  a  hot  plate  which  caused  the  crack  to  propagate  a  distance  through  the  sample.  Only 
samples  with  straight  cracks  that  were  normal  to  the  surface  and  contained  internally  to 
the  sample  dimensions  were  chosen  for  measurement. 

The  V(x)  curves  for  the  circular  and  slot  lenses,  calculated  for  various  defocus 
distances  z,  and  for  the  range  0  <  x  <  500  (im  are  displayed  along  with  measured  V(x) 
curves  in  Figs.  5.1  -  5.2  and  5.3  -  5.4  ,  correspondingly. 

The  theoretical  model  in  Chapter  IV  predicts  that  for  x  near  the  origin,  or, 
equivalently,  lens  approximately  centered  above  the  crack,  there  is  one  cylindrical  (line) 
focus  in  the  lens  rod  located  a  distance  fx  above  the  lens  surface.  Actually,  as  shown  in 
Fig.  5.5,  the  crack  acts  as  a  plane  mirror,  and  the  rays  trace  symmetrically  the  same  paths 
as  they  do  in  the  defect-free  case;  however,  source  and  radiation  points  now  lie  on  the 
same  side  of  the  transducer.  This  results  in  a  ring  (and  not  a  line  )  focus,  as  mentioned 
in  Chapter  II  review  of  the  geometric  ray  method.  Thus  some  disagreement  between 
theory  and  experiment  may  be  expected  for  small  x  in  the  case  of  a  circular  lens.  The 
leaky  rays  entering  the  slot  lens  from  defect-free  sample  form  short  arc  foci  in  the  lens 
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Fig.  5.1  V(  x,  z  =  -300  p.m)  for  a  circular  aperture  lens  :  observed  (sob 
calculated  (dotted)  (scaled  to  match  at  x  =  -500  |im). 
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Fig.  5.2  V(x,z  =  -400  fxm)  for  a  circular  aperture  lens  :  observed  (solid)  and 
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Fig.  5.3  V(x,  z  = -346  Jim)  for  a  slot  aperture  lens  :  observed  (solid)  and  calculated 
(dotted). 
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Fig.  5.4  V(  x,  z  =  -474  Jim)  for  a  slot  aperture  lens :  observed  (solid)  and  calculated 
(dotted). 


rod;  this  error  is  not  expected  to  be  as  significant  as  in  the  case  of  the  circular  lens  since 
for  a  narrow  slot  these  arcs  are  short  and  hence  close  in  shape  to  line  foci. 

When  the  crack  is  near  the  edge  of  the  Rayleigh  circle,  i.e.  x  «  c  =  Izl  tan  0R, 
Rayleigh  waves  are  excited,  scattered  from  the  crack,  and  re-radiated  within  a  small 
distance  and  therefore  the  expression  (4.3.2)  for  coefficient  E  is  no  longer  correct,  as 
mentioned  in  Chapter  IV.  Also,  some  disagreement  is  expected  due  to  the  fact  that  the 
variation  in  the  Rc  and  Tc  with  the  angle  of  incidence  was  ignored. 

5.2  Extracting  the  surface  wave  reflection  coefficient  of  a  crack  from 
V(x)  measured  using  a  slot  lens 

As  mentioned  above,  when  the  lens  is  negatively  defocused  and  the  crack  lies  at 
the  lens  axis,  i.e.  x  =  0,  z  <  0,  the  leaky  rays  trace  paths  that  are  mirror  images  of  leaky 
wave  paths  for  the  defect-free  case.  Therefore,  the  microscope  response  at  that  point  may 
be  represented  in  terms  of  unperturbed  contributions  VG  and  Vl  as 

/VL(z)  VL(z)\ 

V(  x=0,  z)  =  VG(z)  +  Tc  •  VL(z)  +  Rc  •  +  2  (5-2.1) 

where  the  reflected  contribution  has  been  separated  into  two  identical  terms  due  to 
surface  waves  incident  on  the  crack  from  the  right  and  from  the  left.  Note  when  the 
discontinuity  is  not  present,  Tc  =  1  and  Rc  =  0,  and  (5.2.1)  reduces  to  the  usual  form  V(z) 
=  VG(z)  +  VL(z).  At  this  point  fy  tends  to  infinity  (4.3.4),  and  there  is  only  one  focus  in 
the  lens  rod,  located  at  fx  above  the  aperture. 

When  the  lens  is  moved  from  the  crack  in  the  positive  x  direction  (  see  Fig.  4.7), 
the  surface  waves  incident  from  the  left  have  to  propagate  an  additional  distance  of  2x  to 
the  crack  and  back,  thus  gaining  i2kRx  in  phase.  Similarly,  the  rays  incident  from  the 
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right  encounter  the  crack  sooner  and  thereby  lose  i2kRx  in  phase.  Neglecting  other 
variations,  which  will  be  considered  later,  expression  (5.2.1)  can  therefore  be  generalized 
for  x  o  0  as 

Vl(z)  exp(  i2kRx)  +  VL(z)exp(  -i2kRx) 

V(  x,  z)  =  VG(z)  +  Tc  •  VL(z)  +  Rc  • - - 2 - 

=  VG(z)  +  Tc  •  VL(z)  +  Rc  •  VL(z)  cos  2  kRx  .  (5.2.2) 

The  expression  (5.2.2)  can  be  used  as  a  basis  for  the  extraction  of  Rc  from  V(x) 
data  but  first  the  region  of  its  validity  should  be  determined.  Since  fy  varies  as  1/x  for 
small  x,  it  rapidly  approaches  the  plane  of  the  transducer  as  x  increases  from  0.  A  critical 
point  is  encountered  when  this  focal  point  crosses  the  transducer  plane,  which  occurs  at  x 
given  by  (4.3.9).  Evaluating  (4.3.9)  for  the  lens  used  in  this  work  gives  the  location  of 
the  critical  point  as  x  =  15  Jim.  The  amplitude  of  V(x)  calculated  for  z  =  -346  |J.m  is 
plotted  in  Fig.  5.7.  It  can  be  observed  from  the  graph  that  for  0  <  x  <  10  |xm  the  V(x) 
varies  as  a  cosine  function,  in  agreement  with  (5.2.2).  As  discussed  above,  the  model  for 
V(x)  described  in  Chapter  IV  fails  at  x  =  0.  At  that  point,  however,  the  microscope 
response  is  reasonably  well  described  by  (5.2.1).  Thus,  the  expression  (5.2.2)  will  be 
assumed  to  be  valid  for  small  x  and,  therefore,  be  used  for  extracting  a  value  for  Rc. 

The  acoustic  microscope  employed  in  this  work  measures  only  the  amplitude  of 
the  transducer  voltage  IV(x,z)l.  Taking  the  amplitude  of  (5.2.2)  and  squaring  it  for 
convenience  gives 

IIV(x,z)l2  =  IVG  +  Tc  •  VLI2  +  2  IVG  +  Tc  •  VlM  Vl  Rj  •  cos(  2kRx  +  4>(z))  + 

+  IVl  Rj2  cos2  2kRx  .  (5.2.3) 

where  4>(z)  =  arg(VG(z))  -  arg(VL(z))  -  arg  (Rc).  Since  the  measurement  is  carried  out  on 
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a  system  where  the  smallest  translational  step  is  2  Jim,  there  is  a  positional  uncertainty  of 
under  1  |im  in  locating  the  true  x  =  0  point.  The  resultant  error  in  x  can  be  included  as  a 
constant  phase  term  in  <|>  in  (5.2.3). 

When  the  V(x)  measurement  is  performed  at  z  corresponding  to  a  V(z)  maximum 
or  minimum,  the  VG  and  Vl  are  exactly  in  phase  or  180°  out  of  phase,  respectively, 
(V(z)  =  IVGI  ±  IVLI)  and  the  magnitude  of  IVLI  can  be  deduced  by  subtracting  IVGI  from 
V(z).  To  find  IVGI,  the  V(z)  of  teflon  is  scaled  to  match  the  V(z)  of  glass  at  large  z, 
where  IVJ«IVGI. 

The  last  term  in  (5.2.3)  can  be  expanded  as  : 


I  Vl  RJ2  cos2  2kRx 


IVtRJ2 

— 2 —  (1+cos  4  kRx) . 


(5.2.4) 


Evaluating  the  IVLI  and  IVGI  at  several  local  V(z)  maxima  and  minima  gives  IVGI  >  10 
IVLI  in  which  case  the  term  (5.2.4)  is  small  compared  to  the  other  terms  in  (5.2.3)  and 
can  be  neglected.  To  remove  the  x-independent  term  IVG  +  Tc  VLI2  from  (5.2.3),  IV(  0, 
z)l2  -  IV(  x,  z)l2  is  evaluated  as 

IV(  0,  z)l2  -  IV(  x,  z)l2  =  2  l(VG  +  Tc  Vl)  Vl  R<.  I  cos  <J>  (1  -  cos  2kRx) 

-  2  l(VG  +  Tc  VL)  VL  Rg  I  sin  (J)  sin2kRx.  (5.2.5) 


The  objective  here  is  to  find  the  coefficients  multiplying  the  (1  -  cos  2kRx)  and  sin  2kRx 
terms.  When  (5.2.5)  is  evaluated  at  two  points  X!  and  X2,  a  system  of  two  equations  and 
containing  the  two  undetermined  coefficients  A  and  B  results  that  satisfy 


A(1  -  cos  2kRXj)  -  B  sin  2kRxj  =  IV(0,z)l2-IV(  x^,z)l2  and 

A(1  -  cos  2kRx2)  -  B  sin  2kRx2  =  IV(0,z)l2-IV(  x2,z)l2  .  (5.2.6) 
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where  A  and  B  are  independent  of  x  and  given  by 


A  =  2  l(VG  +  Tc  VL)  VL  Rc  I  cos  <}>  and 

B  =  2  l(VG  +  Tc  VL)  VL  Rc  I  sin  (5.2.7) . 

The  system  (5.2.6)  can  be  solved  for  A  and  B  from  the  measured  data.  From  (5.2.7)  it 
follows  that 


VA2  +  B2  =  2I(Vg+TcVl)VlRcI 
and 

<J>  =  tan-1  B/A.  (5.2.8) 

The  factor  l(VG+Tc  VL)I  in  (5.2.7)  is  found  by  noting  that  it  is  exactly  the  square 
root  of  the  x-independent  term  in  (5.2.3)  and  can  therefore  be  obtained  by  substituting 
(5.2.7)  and  (5.2.8)  into  (5.2.3)  for  any  small  x.  Finally,  IRCI  is  found  from  (5.2.8)  using 
the  value  of  IVJ  deduced  from  the  unperturbed  V(z)'s  of  glass  and  teflon,  in  the 
procedure  described  above. 

Since  the  system  of  equations  (5.2.6)  requires  the  knowledge  of  V(x,z)  at  only 
two  x  locations  beside  the  value  V(0,z),  it  is  possible  to  use  the  V(x,z)  data  at  other 
points  to  create  additional  (redundant)  equations  and  simultaneously  solve  any  two  of 
them.  The  resulting  values  of  IRJ  can  then  be  checked  for  consistency  and  if  in 
reasonable  agreement  can  then  be  averaged. 

The  results  of  applying  this  algorithm  to  extract  IRCI  from  V(x)  measured  at 
different  defocus  distances  z  are  presented  in  Table  I.  Since  the  signal  detected  by  the 
transducer  weakens  significantly  as  Izl  increases,  only  V(x,z)  measured  at  small  Izl  are 
deemed  reliable,  i.e.,  least  corrupted  by  noise.  It  is  seen  that  the  values  of  IRCI  in  Table  I 
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for  z  >  -210  jxm  compare  favorably  with  the  value  IRCI  =  0.54  deduced  by  Achenbach,  et. 
al.  [14]  for  a  deep  crack.  Furthermore,  the  average  value  of  IRCI,  calculated  for  all  z  in 
Table  I,  is  found  to  be  0.55. 

Angel  and  Achenbach  [3]  have  also  produced  a  study  of  the  dependence  of  IRCI 
on  the  crack  depth,  reproduced  in  Fig.  4.16.  Such  a  correlation  could  be  used  to 
determine  the  crack  depth  once  IRCI  is  found.  Thus,  the  V(x,z)  can  be  used  as  an 
important  tool  in  defect  characterization. 
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z(|xm) 

Averaged  IRCI 

-72 

0.50 

-130 

0.53 

-208 

0.56 

-346 

0.60 

-416 

0.75 

-474 

0.34 

-534 

0.66 

-604 

0.48 

Table  I.  Values  of  crack  reflection  coefficient  IRCI  deduced  from  V (x)  measured  at 
various  defocus  distances  z. 
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Appendix  A 


Derivation  of  geometric  quantities  pR,  fx,  fy,  pRmin»  pRmax>  cmin 111(1  cmax- 

A.1  Derivation  of  pR 

The  leaky  ray  that  appears  to  originate  from  the  lens  focal  point  crosses  the  lens 
surface  at  a  distance  p  =  Pr  from  the  lens  axis.  It  follows  from  Fig.  A.l  that: 

Pr  =  -pR^)tan0R.  (A.l) 

Solving  for  pR  gives: 

Pr  =  (f  -  R)  sin  0r  cos  0r  +  sin  0r  (f  -  r)2  cos2  0R  -  (f2  -  2Rf)  .  (A.2) 

Substituting  0R  =  28°  (glass),  R  =  3  mm  and  f  =  4mm  into  (A.2)  gives 
Pr=  1.8  mm. 

A.2  Derivation  of  fx 

The  leaky  rays  strike  the  lens  as  parallel  family  of  rays  and  focus  at  a  distance  fx 
above  the  lens  surface.  In  Fig.  A.2  a  leaky  ray  through  the  lens  center  of  curvature  and 
the  principal  ray,  that  appears  to  come  from  the  lens  focus  are  shown.  It  follows  from 
this  Figure  that: 
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f*=dre£-VR2-pR2  • 


(A.  3) 


A.3  Derivation  of  fy 

The  leaky  rays  re-radiated  from  the  arc  BAC  (Fig.  4.8)  form  a  cylindrical  (line) 
focus  at  a  distance  fy  above  the  lens  surface.  These  rays  appear  to  originate  from  a  point 
shown  in  Fig.  4.9  re-drawn  for  convenience  in  Fig.  A.3.  This  line  focus  occurs  at  the 
intersection  between  the  principal  ray  and  a  ray  drawn  through  the  lens  center  of 
curvature.  Using  Fig.  A.3  ,  fy  can  be  found  from 


2x  f  -  R  +  2x/tan  0R 
Pr  "  fy+^R2-pR2  ' 


(A.4) 


Solving  for  fy: 


fy  = 


(f-R)  Pr 
2x 


+  f. 


x » 


where  fx  is  given  by  (A.3). 


(A.5) 
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A.  4  Derivation  of  pRmin,  pRinax,  cmin»  cmax 


The  limiting  rays  of  the  leaky  ray  bundle  that  strikes  the  transducer  cross  the  lens 
surface  at  PRmin  and  PRma*,  shown  in  Fig.  A.4.  Using  fx  «  D,  these  quantities  can  be 
approximated  as: 


pRmin  =  PR-(RT-pR)fx/°  * 

PRmax  =  PR+ (RT  +  Pr)^/0  •  (A-6) 

These  limiting  rays  appear  to  originate  from  points  at  p  =  cmjn  and  p  =  cmax, 
respectively,  on  the  sample  surface, as  shown  in  Fig.  A.4.  In  Fig.  A.4  the  lens  surface  in 
the  vicinity  of  p  =  Pr  is  approximated  as  a  straight  line  at  an  angle  a  with  respect  to  a 
vertical  line.  The  angle  a  is  found  from 

a  =  cos*1  ( pR  /  R ) .  (A.7) 

The  quantities  lzi„inl  and  l^^l,  shown  in  Fig.  A.4  are  found  to  be 


IZmin1  =  lzl  -  (PR-pRmin)  (cot  0r  +  cot  a) 

'zmaJ  =  lzl  +  (pRmax'PR)  (cot  0R  +  cot  a)  •  (A-8) 

It  follows  from  Fig.  A.4  that 

Cmin  =  lzmin!  tan  0R 

cmax  =  ^Zmax^  tan  0R .  (A.9) 
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